Topological versus rheological entanglement length in primitive path analysis protocols 
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Primitive path analysis algorithms are now routinely employed to analyze entanglements in com- 
puter simulations of polymeric systems, but different analysis protocols result in different estimates 
of the entanglement length, A^'^. Here we argue that standard PPA measures the rheological entan- 
glement length, Nl^'^° , typically employed by tube models and relevant to quantitative comparisons 
with experiment, while codes like Z or CRcTA also determine the topological entanglement length, 
NI°^° . For loosely entangled systems, a simple analogy between between phantom networks and the 
mesh of entangled primitive paths suggests N^^''° « 2Nl°^° . This result is in excellent agreement 
with reported values for poly-ethylene, poly-butadiene and bead-spring polymer melts. 

PACS numbers: 83.80.Sg (Polymer melts), 83.10.Rs (MD simulation rheology), 61. 25. he (liquid structure 
polymer melt) 
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Understanding the role of entanglements in the dy- 
namics of high molecular weight polymeric liquids is a 
classical subject of polymer physics [ll, Q • Tube models 
[ll-t3)i describe the neighborhood of a single chain in terms 
of a mean-field confinement potential, which transiently 
localizes thermally fiuctuating chains in a tube-like re- 
gion in space. Primitive paths were originally introduced 
in a thought experiment to determine the tubes confin- 
ing individual polymers in an entangled polymer melt 
or network. The idea is to identify the random walk- 
like tube axis with the shortest ("primitive") paths be- 
tween the end points of the original chains into which the 
chain contours can be contracted without crossing each 
other. Corresponding primitive path analysis (PPA) al- 
gorithms |4l4lC| have become a standard tool in computer 
simulations of polymeric systems. In particular, we have 
claimed that the primitive path analysis (i) endows the 
tube models with predictive power for dynamical pro- 
cesses in flexible chain melts [ll| and (ii) provides a quan- 
titative unified picture of chain entanglement over a much 
wider class of polymeric systems jTo|. However, estimates 
of the melt entanglement lengths based on different PPA 
variants seem to vary within a factor of two . Is there 
something amiss? 

The present note was triggered by a recent study by 
Hoy et al. [l^l , which contains a comparison of PPA and 
Z results for polyethylene and bead-spring polymer melts. 
The original PPA algorithm 0, Q replaces the intra- 
chain interactions by harmonic bonds between neighbor- 
ing beads. In the course of an energy minimization with 
fixed chains ends, the bond springs try to reduce the bond 
length to zero and pull the chains taut. The inter-chain 
excluded volume interactions ensure that different chains 
do not cut through each other. In contrast, Kroger's 
numerically much more efficient Z-code uses geometrical 
operations to shorten primitive paths, which are treated 
as tensionless lines rather than multibead chains. Other 
PPA algorithms 

follow similar recipes. 

A first observation by Hoy et al. [l2| well worth repeat- 
ing is that there is no substantial difference between the 



results from the PPA and the Z algorithms, if one com- 
pares identical observables (typically the primitive path 

length, {Lpp) or ^ {Lpp})- Nevertheless, their results sug- 
gest a relation of N^^^ « 2A^^ between entanglements 
lengths inferred following the PPA and Z protocols: PPA 
practitioners measure the number of monomers per PP 
Kuhn length, Opp, while is the number of monomers 
between interior kinks at positions, where two primi- 
tive paths intersect. That < N^^^ and that there 
is hence an orientational correlation between primitive 
path segments was first pointed out by Tzoumanekas and 
Theodorou 0. For loosely entaiigled systems, the results 
should agree on a scaling level [l3|- But why is there a 
factor of two? And are we sure, that it is really N^^^ 
and not which should be compared to experimentally 
measured entanglement lengths? 

In the following, we will argue that these questions are 
easily resolved by considering a corresponding problem 
in the framework of the classical theory of rubber elas- 
ticity [l^, which deals with the behavior of networks of 
Gaussian chains under strain. In particular, we will find 
that the present issue is related to the "frontfactor" de- 
bate |14| - [l7| for the predicted elastic moduli. Theiunc- 
tion affine model [la| and the phantom model [ij, [l^ 
both predict shear moduli of G = jPstranciksT, where 
Pstrand = p/Ns dcuotc the number density of elastically 
active network strands and monomers respectively. But 
while 7=1 for the junction affine model, the phantom 
model predicts a dependence 7 = (1 — 2//) on the func- 
tionality of the junction points (or, more appropriately, 
the cycle rank J7|] of the network). It took until 1976 for 
Flory [17| to resolve the discussion. The salient points 
are apparent from the following simplified argument : 
in a phantom network, the thermal fiuctuation of a junc- 
tion point around its rest position are of the order of 
Ngb'^/f, where (i?^) — Ngb^ is the mean-square end- 
to-end distance of the network strands of length Ng in 
a corresponding linear chain melt. If one assumes that 
neighboring junction points fluctuate independently (ob- 
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viously a simplification) , then the thermal fluctuations of 
the extension of a network strand are given by 2/ f{E?). 
Since instantaneous random crosshnking does not change 
the static ensemble average of the strand end-to-end dis- 
tances, one can thus deduce the quenched part of the 
fluctuations: the mean-square average strand extension 
in the network is (1 - 2//)(i?2) = 7(i?2). The same 
prefactor appears in the shear modulus predicted by the 
phantom model, 
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(1) 



because the thermal fluctuations of a phantom network 
are invariant under strain, while the average strand 
extensions undergo an affine transformation with the 
macroscopic deformation The earlier junction affine 
model neglects this subtle mathematical point and de- 
rives the elastic reponse from an affine transformation of 
the distribution of instantaneous end-to-end distances, 
so that 7=1. In real networks, the situation is inter- 
mediate between these two limits, since entanglements 
contribute to the junction localization [2T|-[23j. 

For a phantom network prevented from collapse by pe- 
riodic boundary conditions or surface anchored junction 
points, PPA converges to the ground state with all junc- 
tion points localized at their average positions [l^. To 
understand the relation between the above consideration 
and the interpretation of PPA results, we need to ana- 
lyze the statistics of (primitive) paths through the net- 
work (these paths could, for example, correspond to long, 
multiply crosslinked precursor chains |25i, ^] ) . On large 
scales, N ^ Ng, the paths follow the same random walk 
statistics as the original chains with (i?^)(iV) = Nb'^. 
In contrast, PPA smoothens out local thermal fluctua- 
tions on the mesh scale, so that the paths become locally 
straight between junction points. By construction, the 



Kuhn length. 



of the paths is given by 



bN, 



Lpp{N) ^(1 - 2//)A^, 



(2) 



where Lr, 



-- ^^/{l - 2/f){R^){Ns) is the contour 
length distance along the path and where we have 
deduced the mean-square contour length per network 
strand from the above arguments for phantom networks. 
In analogy to the standard PPA protocol, we now con- 
sider the number, N^^^ , of monomers per Kuhn length 
of the primitive path: 
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^pp 



Lpp/N 1 - 2/f 



(3) 



In particular, we see by comparison with Eq. ([T]) that 
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Neglecting the subtle difference between slip-links [27[ 
and crosslinks, we see that the number of monomers be- 
tween entanglements or topological entanglement length, 
]\[topo _ j\j-z direct analogue of the strand length 

Ns in a polymer network. Nevertheless, it is the less in- 
tuitive rheological entanglement length, N^^'^" = N^^^, 
which should be compared to experimentally determined 
values of , because the latter are derived from the 
melt plateau modulus (or the closely related entangle- 
ment modulus 0) and the relation Ge = pkBT/Ng. 

The reported values of N^''^° /N^''p° = 2.6 
and N^'"'°/Nl°P° = 1.9 and 2.0 [H for poly- 
ethylene, Nl'^<'°/Nl°P° = 2.1 for poly-butadiene and 
jyrheo/jytopo ^ j^ g ^ f^^. Krcmer-Grest melts [H are 

compatible with the interpretation of entanglements as 
binary contacts (or / = 4 functional sliplinks) between 
primitive paths. We emphasize, that the present anal- 
ogy only holds for fully flexible, loosely entangled chains. 
Tightly entangled chains fall into a different regime with 

jSjrheo ^ j^topo 
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